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ABSTRACT 

Context. The expansion of Galactic H ii regions can trigger the formation of a new generation of stars. However, little is know about 
the physical conditions that prevail in these regions. 

Aims. We study the physical conditions that prevail in specific zones towards expanding H ii regions that trace representative media 
such as the photodissociation region, the ionized region, and condensations with and without ongoing star formation. 

Methods. We use the SPIRE Fourier Transform Spectrometer (FTS) on board Herschel to observe the Hii region RCW 120. 
Continuum and lines are observed in the 190 - 670 pm range. Line intensities and line ratios are obtained and used as physical 
diagnostics of the gas. We used the Meudon PDR code and the RADEX code to derive the gas density and the radiation field at nine 
distinct positions including the PDR surface and regions with and without star-formation activity. 

Results. For the different regions we detect the atomic lines [Nil] at 205 pm and [Cl] at 370 and 609 pm, the *^CO ladder between 
the y = 4 and 7=13 levels and the *^CO ladder between the 7 = 5 and 7 = 14 levels, as well as CH^ in absorption. We find gas 
temperatures in the range 45 - 250 K for densities of 10“* - 10® cm“^, and a high column density on the order of Ah ~ 10^^ cm“^ 
that is in agreement with dust analysis. The ubiquitousness of the atomic and CH^ emission suggests the presence of a low-density 
PDR throughout RCW 120. High-excitation lines of CO indicate the presence of irradiated dense structures or small dense clumps 
containing young stellar objects, while we also find a less dense medium (Ah ~ 10^** cm“^) with high temperatures (80 - 200 K). 

Key words. ISM: Hii regions - ISM: individual objects: RCW 120 - ISM: bubbles - ISM: photon-dominated region (PDR) - 
infrared: ISM - submillimeter: ISM 


1. Introduction 

Star formation occurs on the borders of Galactic Hii regions. 
Different physical processes may be at work to trigger this star 
formation ([Elmegreen & Lada |1 977 [ Deharveng et al. 2010|l. 


This phenomenon has been studied in detail over the past ten 


years 

([Deharveng et al.[[2003[ 

Zavagno et al. 

2007 Pomares 

et al. 

2009| Zavagno et al.|2010 

Paron et al.|20] 

1 

Davies et al. 

2012 

1. All these studies concentrate on specific Ga 

actic HII re- 


gions where young stellar objects (YSOs) are observed on their 
borders. This phenomenon is also observed in nearby galaxies 
(see, e.g., the case of Nil in the Large Magellanic Cloud). The 
Spitzer satellite and the GLIMPSE and MIPSGAL surveys of 
the Galactic Plane have revealed that we are living in a bubbling 
galactic disk where H ii regions have a clear impact on their envi¬ 
ronment. A study using Spitzer GLIMPSE and MIPSGAL data 
combined with ATLASGAL data on 102 bubbles have shown 
that the star formation triggered by H ii regions is an important 
phenomenon in our Galaxy. Up to 25% of the bubbles show trig¬ 
gered massive-star formation on their border (|Deharveng et al. 

[2OTOI 1. 

The Herschel satellite offers a unique opportunity to study 
the star formation triggered by Galactic Hit regions. Thanks to 
its sensitivity and its large wavelength coverage in the far in¬ 


* Herschel is an ESA space observatory with science instruments 

provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 


frared, Herschel is perfectly suited to the study of the earliest 
phases of star formation. 

We are engaged in several guarante ed (HOBYS,|Motte et al. 
|2010[ “Evolution of Interstellar Dust”, [Abergel et al.|2010[ ) and 
open time (Hi-GAL, [Molinari et al. 20101 key programs on 
Herschel that aim to characterize this way of forming stars. We 
used PACS and SPIRE photometers to characterize the emission 
in the 60 - 500yum range towards Galactic Hit regions. This al¬ 
lows us to detect and characterize the properties of young stellar 
objects (YSOs) observed towards these regions. We also used 
the PACS and SPIRE spectrometers to derive the physical con¬ 
ditions towards these regions. 

Here we present results of Herschel SPIRE-ETS spec¬ 
troscopy obtained towards the Galactic Hii region RCW 120. 
This region is among the closest H ii regions with a distance of 
~ 1.3 kpc ( |Russeil|2003 1, which combined to its angular diame¬ 
ter of ~ 7'.5 Anderson et al. ( 2015| l results in a physical size of 
~ 3pc. RCW 120 is ionized by the single star CD - 38°11636, 
with a spectral type 06-8V/III according to the latest measure¬ 
ments by jMartins et al.| ( [2010| l. 

Labeled “the perfect bubble” throughout the literature, recent 
single-dish observations of the lowest CO transitions (Anderson 
et al.||2015| Torii et aT||2015| l fail to detect an expanding shell 
of fore- and background material, which would indicate a 3-D 
structure. However, [Anderson et al.| ( |2015] l finds discrete “holes” 
in the PDR, through which the ionizing radiation is escaping, 
and [Torii et ah] ([2015|l propose an explanation of the observed 
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Fig. 1. Herschel SPIRE-FTS pointing towards RCW 120; Herschel SPIRE 350 image (background) of RCW 120 on which the 
SEW (blue) and SSW detectors (yellow) are superimposed. The nine pointings observed are labeled according to their location (see 
text). 


morphology of RCW 120 through the collision of two clouds, 
which formed the ionizing star. 


by the collect and collapse process on the border of an H n re¬ 
gion. 


Observations in mm and sub-mm wavelengths ( [Zavagno] 
et al. 2007 [ Deharveng et al.| 2009|l show a fragmented neu¬ 


tral layer along its PDR, which contains five of the eight mm- 
condensations found by [Zavagno et al.| ( |2007| l. According to 
Deharveng et al. ( 2009[l the mass of the layer is ~ 2000 Mq. 


The most massive condensation has been resolved into a chain 


of several Class I objects dDeharveng et al.|2009 

1 , likely an ex- 

ample of Jeans instability, while Zavagno et al. ( 

20101 found a 


massive 8 - 10 Mq YSO towards this condensation, which they 
suggest is the first detection of a massive Class 0 object formed 


2. Observations and data reduction 


RCW 120 was observed in spectroscopy with the SPIRE-Fourier 
Transform Spectrometer (SPIRE-FTS, Griffin et al.|2010[ ) on 26 
February 2010 (O P 288), as part of the Evolution of Interstellar 
Dust key program (Abergel et al.[2010']p' The aim of the program 


* The reduced data cubes are available on the HESIOD portal (http: 
//idoc-herschel.ias.u-psud.fr/sitools/client-user/I 


2 











































Rodon, J.A. et al.: Herschel SPIRE-FTS observations of RCW 120 


Table 1. SPIRE-FTS observations on Operational Day 288 (February 26 2010). 


Pointing 

RA (J2000) 

Dec (J2000) 

Obs. ID 

Start time (UT) 

Int. time (s) 

Repetitions 

C6 

17:11:47.89 

-38:19:48.2 

1342191225 

15:21:40.00 

1353 

8 

CLay 

17:12:04.91 

-38:24:03.5 

1342191226 

15:44:38.00 

1353 

8 

Cl 

17:12:08.93 

-38:30:43.1 

1342191227 

16:07:38.00 

543 

2 

C8 

17:12:19.81 

-38:34:04.0 

1342191228 

16:17:04.00 

543 

2 

IRS2 

17:12:20.77 

-38:18:21.1 

1342191229 

16:26:43.00 

1353 

8 

HII 

17:12:24.47 

-38:27:14.3 

1342191230 

16:49:45.00 

543 

2 

C2 

17:12:34.15 

-38:30:51.9 

1342191231 

16:59:11.00 

543 

2 

IRSl 

17:12:40.89 

-38:27:07.9 

1342191232 

17:08:37.00 

543 

2 

OFF 

17:13:01.35 

-38:27:13.3 

1342191233 

17:18:03.00 

1353 

8 


Notes. Coordinates are for the central pixels of the SLW and SSW arrays (pixels C3 and D4, repectively). 


is to study the physical conditions that exist in regions of trig¬ 
gered star formation. For this purpose we selected 8 representa¬ 
tive positions towards RCW 120: towards the ionized gas (point¬ 
ing HII), towards the photo-dissociation region (PDR) without 
any condensation or star formation (pointing CLay), and towards 
condensations with ongoing star formation (see Sec. 2.11. In 
RCW 120, nine condensations suspected of harboring ongoing 


star formation were observed in the mm continuum by Zavagno 
|et al.| ( |2007| l. For the FTS observations, we selected six such 
condensations, labeled Cl, C2, C6, C8, IRSl and IRS2. A posi¬ 
tion “off” the Hit region was also observed (pointing OFF). All 
nine pointings are shown in Fig.[2 The short wavelength (SSW, 
194-320 yum) and long wavelength (SLW, 313-671 /rm) amys 
of the SPIRE-FTS are overlaid on the 350 fim SPIRE image of 
RCW 120, in yellow and blue, respectively. The arrays are com¬ 
posed of an hexagonal grid of 19 pixels for SLW, and 37 pixels 
for SSW. The pixel distribution and labels are shown in Fig. 
There are two “dead” pixels in SSW, pixels F4 and D5, and are 
thus not shown in Fig. T](missing yellow circles). 

All the observations were performed in the high-resolution 
mode of the SPIRE-FTS. In general two scan repetitions were 
done, except for positions expected to be fainter, where eight 
scan repetitions were taken. An observation summary is given in 
Tabled 
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Fig. 2. Disposition of the detectors (pixels) on the SSW (right) 
and SLW (left) arrays of the SPIRE-FTS instrument for our ob¬ 
servations as seen projected in the sky. The two gray pixels in 
the SSW array mark dead pixels. 


The data reduction was performed using the Herschel 
Interactive Processing Environment (HIPE) version 9.0 build 
3048 with the calibration tree 9_1. An iterative spectral line fit¬ 
ting routine has been used to extract the line parameters (line 


center, intensity and associated errors) from the unapodized FTS 
spectra. The routine used the unapodized instrumental line shape 
(the classical sine function) as the lines were checked to be un¬ 
resolved by the FTS. The two scan directions were treated sep¬ 
arately. As a conservative approach, the retained final error is 
the maximum value between the fitting procedure errors and the 
uncertainties obtained from the two direction results. 


2.1. Pointings 


The nine pointings were selected on the basis of what was known 


from the infrared and millimeter study of RCW 120 by Zavagno 
|et al.| ( |2007| l. Their 1.3 mm continuum emission map revealed 
the presence of nine condensations, five of which are located on 
the PDR surrounding the ionized region of RCW 120. The idea 
of this SPIRE-FTS study was to study the physical conditions 
that prevail in the different regions located around RCW 120, 
regions that sample different evolution stage of star formation or 
different typical media (ionized region, PDR). 

One FTS pointing was obtained on each of the main conden¬ 
sations (Cl, C2, C6, C8 and IRSl). The stellar content of these 


condensations is described in detail in Zavagno et al. (20071 and 
Deharveng et al.|(2009|l. The properties of the point sources 


observed in this region with Herschel PACS and SPIRE are de¬ 
scribed in Zavagno et al. ( 2010| l. An off position (OFF) was ob¬ 
tained to serve as a reference. A position was obtained towards 
the ionized gas at the center of the ionized region (HII). A point¬ 
ing was obtained towards the PDR, on the northeastern side of 
the region (CLay). All these pointings are described in more de¬ 
tail below. 

Cl: the most dense and massive condensation observed on 
the borders of RCW 120. This region contains a massive Class 0 
source revealed with Hersc hel ([Zavagno et al. [2070) 1 and a chain 
of young, lower mass stars ( Deharveng et al.|2009| l. 

C2: this condensation contains a bright Class I source and is, 
as Cl, also located on the southern (and densest) border of the 
PDR. 

C6: this condensation is located in the northeastern part of 
the region but not in direct contact with the main ionizing front. 
However, a deep Ha image shows that the ionized gas leaking 
from the H ii region hits this zone and is in direct contact with it 
(see Fig. 3 in Deharveng et al.|2009| l, possibly acting on the star 
formation there. This condensation is dense and contains many 
bright Class I YSOs. 

C8: this condensation is located in the southern part of the 
ionized region and has probably been shaped by the leaking UV 


3 










































Rodon, J.A. et al.: Herschel SPIRE-FTS observations of RCW 120 


■$ 

N 

■55 

E 



300 400 500 600 


Wavelength (jim) 


Wavenumber (cm’’) 

30 25 20 15 



Fig. 3. SPIRE-FTS SLW spectra towards the richest pointings, 
at the position of the central pixels (see Sec.[^. The spectra are 
normalized at the brightest line on each spectrum. The lines de¬ 
tected are marked and labeled. 


field (see Fig. 16 in |Deharveng et al.|2009| l. This region contains 
YSOs seen by Herschel and is a site of active star formation. 

IRS 1: this region is located in the west part of the PDR and is 
labeled condensation 4 in Zavagno et al. ( 2007| l. It was renamed 
IRSl for the FTS pointing because this region hosts bright IR 
YSOs observed with Spitzer, but the peak of the 1.3 mm map 
was devoid of sources (see Fig. 8 in Zavagno et al. [2007) 1. This 
region is of interest because it samples very different physical 
conditions on a small spatial scale. 

IRS2: this region is located on the northern part of the PDR. 
The central pixel position was chosen to match with the position 
of the source IRAS 17089 - 3814. 


CLay: with the idea of sampling all the different physical 
conditions that exist towards RCW 120, we selected this zone 
known to be free of condensations detected in the mm map and 
star formation. Therefore this region should be representative of 
the PDR. At the time of the selection, the only process envi¬ 


sioned to form the layer was the collect and collapse (Elmegreen 
& Lada|1977|l, therefore the name Collected Layer (CLay). 


HII: this pointing is obtained towards the ionized region. No 
star formation and no 1.3 mm continuum emission is observed 
towards this region. 

OFF: the off position was chosen by looking at the Spitzer 
mid-IR map of RCW 120. The PACS and SPIRE maps confirm 
that this is a low-emission region. However, after the analysis we 
find that this position has similar emission to other parts of the 
region (e.g., CLay, see Sec.[^. 


Fig. 4. Same as Fig. [^but for the poorest pointings. 


Wavenumber (cm’’) 

50 45 40 35 



Fig. 5. Same as Fig.j^but for all SPIRE-FTS SSW pointings. 
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Fig. 6. SPIRE-FTS SLW and SS W combined spectra for pointings Cl (left) and IRS 1 (right), showing how in some cases the overlap 
between detector arrays is not perfect (see text). 


3. Results 


3.1. Measured spectra 


Each of the pixels of SLW and SSW arrays produces one spec¬ 
trum, therefore we obtained 486 spectra. Owing to space con¬ 
straints, we will only show a few selected spectra in this paper, 
in Figs. 00 and These are the spectra at the central pixels 
of SLW and SSW for each pointing, which are the only cases 
where a spectrum on the whole wavelength range of the FTS at 
a given pointing can be obtained since the central pixels are the 
only ones whose centers are spatially coincident on the sky. For 
other pixels there certainly is a spatial overlap between the SLW 
and SSW (see Figs. and [^, however the relative centers are 
shifted. 

The lines detected are marked in Figs. toThese are the 
main lines observed in the obtained spectra, and the detections 
are discussed in the next section. 


The overlap between the SLW and SSW portions of the spec¬ 
tra is not perfect. It can be seen in Fig. that in some cases 
there is a vertical offset between the two portions of the spectra. 
This offset is due to the morphology of the regions mapped, and 
the complex properties of the SPI RE beam (see, e.g., |Makiwa | 
et al.|2tfr3| Spinoglio et al. 2012[ Fletcher, L. N. et al.| 2012 


and references therein). The flux calibration for the FTS obser¬ 
vations uses a Relative Spectral Response Function based on the 
telescope model emission. By default it is assumed that the ob¬ 
served source is uniformly extended over the entire beam. Any 
source morphology (compact-like) would affect the calibration. 
In our case, the result is that a relatively compact source (e.g.. 
Cl) will produce an apparently broken continuum, while for 
an extended source (e.g., IRSl) the detected signals are sim¬ 
ilar in both detectors. New calibration tools have been devel¬ 
oped for semi-extended sources (see [Etxaluze et al. 2013[ Wu 
|et al.|20l3] l. However, as we do not know the morphology of our 
source (which could vary with frequency), we are not able to use 


these tools. In Section 3.3 we address the consequences this has 
on our results. 

Because the FTS is a Fourier spectrograph, the line profiles 
have the characteristic sine function shape. This is more no¬ 
ticeable in the strong lines, for example the CO(6 - 5) line at 
pointing IRSl, shown in Figs.l^andl ^ The spectral resolutio n is 
0.048 cm ' throughout the FTS ban ^Swinyard et sT. 2010| . 


Wavenumber (cm ') 

28.4 28.2 28.0 27.8 27.6 27.4 



Fig. 7. Zoom of the normalized CH^(1 - 0) absorption feature 
averaged over all pixels for each pointing. 


3.2. Detected line intensities 

In the different spectra we detect the '^CO ladder from the J = 
4 —> 3 to the 7=13^12 transitions; the '^CO ladder between 
the 7 = 5^4 and the 7 = 8 —> 7 transitions; the [Cl] ^Pi Pq 
and ^P 2 Pi transitions ([CI](1 - 0) and [CI](2 - 1) from now 
on); and the [Nil] Pq transition ([NII](1 - 0) from now 

on). In Tables andare the details of the lines detected in 
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Table 2. Properties of the CO lines detected in emission in the central spectrum of each pointing. 



4-3 

5-4 

6-5 

7-6 

CO 

8-7 9-8 

10-9 

11-10 

12-11 

13-12 

A 

(A'm) 

650.3 

520.3 

433.6 

371.7 

325.2 

289.1 

260.2 

236.6 

216.9 

200.3 

F 

‘-'up 

(K) 

55.3 

S3 

116.2 

154.9 

199.1 

248.9 

304.2 

365 

431.3 

503.1 

FWHM (") 

40.4 

32.6 

29.4 

34.8 

36.8 

19.2 

17.7 

17.6 

17 

16.8 

Cl 


2.40 ±0.32 

6.33 ±0.62 

16.04 ±0.47 

25.59 ± 1.77 

37.20 ±0.64 

51.76 ±3.88 

72.53 ±4.30 

77.53 ±4.73 

66.87±5.17 

39.96 ±5.60 

C2 


2.15 ±0.27 

4.24 ±0.23 

8.21 ±0.16 

12.60± 1.38 

17.43 ±0.29 

17.47 ±0.76 

16.56± 1.58 

11.47 ±0.92 

6.08 ± 1.01 

5.90 ±1.29 

C6 


0.73 ± 0.05 

1.15±0.06 

1.36 ±0.05 

1.45 ±0.49 

1.36±0.18 






C8 

's 

0.92 ±0.25 

1.55 ±0.27 

2.18±0.17 

2.02 ±0.42 

1.37 ±0.40 






IRSl 

7 

1.44 ±0.06 

3.23 ±0.39 

6.71 ±0.23 

10.76 ±0.92 

13.36±0.12 

12.92 ±1.74 

11.61 ±0.56 

7.77 ± 1.62 

4.06 ± 1.19 


IRS2 

60 

0.43 ± 0.06 

0.50±0.10 

0.29 ±0.07 








HII 

7 

0.37 ±0.05 










CLay 


0.22 ±0.02 










OFF 


0.53 ± 0.06 

0.51 ±0.04 

0.30 ±0.07 









emission on each of the central spectra, shown in Figs. [^I^and 

HI 

In the pointings towards a condensation (Cl, C2, C6, C8, 
IRSl, and IRS2) or the PDR (CLay) of RCW 120, the number 
of lines detected on each pixel correlate with the position of said 
pixel in the sky in such a way that the pixels towards the PDR 
material or the condensation show richer spectra than the pixels 
“off” the PDR or the condensation. This means that for exam¬ 
ple in pointing Cl pixels Cl, Dl, and El show fewer lines than 
the others. In the other pointings, towards diffuse gas (HII and 
OFF), the lines detected are mostly the same for all pixels of a 
given pointing. The only line that is detected in every pixel of 
every pointing is the [NII](1 - 0) transition. The [Cl] transitions 
are present in almost all the pixels. This is also true for the first 
three CO transitions detected, except in pointings HII and CFay, 
where only the lowest CO transition detected, i.e., 7 = 4^3, 
is present. The ^^CO lines are detected only in pointings Cl, C2, 
and IRS 1. Those pointings are also the only ones with detections 
of the higher CO transitions. 

The CH^(1 - 0) line has been also detected, but in all the 
cases it appears in absorption. This CH^ transition is seen in 
emission in most of the Orion Bar ( Naylor et al.|2010 1 while it 
is seen in absorption in two ultracompact H ii (UCH ii) regions 
located near the Galactic plane ( Kirk et al.|2010| l. Therefore we 
can expect that the corresponding lines seen in the direction of 
RCW 120 should be the result of both emission and absorption 
mechanisms. Because of the low spectral resolution of the FTS 
these lines are not spectrally resolved and we cannot obtain their 
equivalent width, thus no quantitative results can be derived from 
the FTS spectra. Nevertheless, we can see from Fig. that the 
largest absorption seems to be in the direction of the OFF po¬ 
sition, while the pointings towards the PDR of RCW 120 show 
the weakest signatures. All this suggests that CH^(1 - 0) would 
be tracing diffuse gas along the line of sight, while RCW 120 
is in fact emitting in this line, dampening the absorption in the 
pointings towards the H ii region (see, e.g.,|Falgarone et al.|2010 
INagy et al.|20T3] l. 


3.3. Determination of physicai properties using RADEX 

In Fig.j^we present the integrated line intensities of the different 
positions for ^^CO (crosses) and '^CO (plus signs). 

Comparing the curves of the integrated line intensities, four 
groups can be distinguished; 


Table 3. Same as Table |2but for '^CO 




5-4 

13 

6-5 

CO 

7-6 

8-7 

A 

(um) 

544.2 

453.5 

388.7 

340.2 

E„p 

(K) 

79.3 

111.1 

148.1 

190.4 

FWHM (") 

32.9 

30 

34 

36.1 

Cl 


1.48 ±0.37 

2.35 ±0.45 

3.74 ±0.52 

3.97 ±0.60 

C2 


0.86 ±0.09 

1.23 ±0.31 

1.54±0.13 

0.94 ±0.31 

C6 

rt 





C8 

B 

O 





IRSl 


0.73 ±0.16 

0.71 ±0.28 

0.78±0.10 

0.65 ±0.11 

IRS2 






HII 

O 





CLay 






OFF 







Table 4. Same as Table|^but for atomic species 



[Cl] 

1-0 2-1 

[Nil] 

1-0 

A (pm) 

609.1 

370.4 

205.2 

E„p (K) 

23.6 

62.5 

70.1 

FWHM (") 

37.2 

34.8 

16.9 

Cl 

L12±0.34 6.80±0.57 

84.99 ±5.46 

C2 ^ 

1.01 ±0.28 5.50 ±0.20 

46.04 ±3.87 

C6 

0.61 ±0.06 1.92 ±0.06 

10.02 ±0.39 

n 

00 

cm' 

0.67 ±0.19 1.66 ±0.11 

15.06 ± 1.02 

IRSl 7 

0.95 ±0.14 3.72 ±0.26 

26.33 ±0.72 

IRS2 £? 

0.53 ±0.06 1.18 ±0.05 

13.27 ± 1.45 

HII 7 

o 

0.48 ±0.05 0.88 ±0.11 

34.39 ± 1.24 

CLay C 

0.29 ±0.03 0.87 ±0.09 

13.13± 1.20 

OFF 

0.53 ±0.07 0.98 ±0.08 

19.96 ±0.34 


1. Cl, C2 and IRSl show high-excited '^CO and '^CO lines 
and a small ^^CO/'^CO ratio for low-excited lines, 

2. C6 and C8 show no '^CO lines, no high-excited '^CO lines 
and the low-excited '^CO lines have a moderately high in¬ 
tensity, 

3. IRS2 and OFF have no ^^CO lines, no high-excited '^CO 
lines and the low-excited '^CO lines have a rather low inten¬ 
sity, and 

4. CFay and HII show only one, the '^CO(4 - 3) line. 
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Fig. 8. Fit of integrated line intensities of the observed '^CO and '^CO lines (symbols) with RADEX for three densities: n - 10"^ 
(green lines), n - 10^ (blue lines) and n - 10® (black lines). The resulting temperatures and CO column densities are shown 
in Table |5] 


In order to asses the kinetic temperature, density and column 
density in which the CO lines arise in the different positions, 
we fit the integrated line intensities of the observed '^CO and 
'^CO lines with the non-LTE and (local) radiative transfer code 
RADEJ^ We use a grid of input parameters, where the gas den¬ 
sity varies between lO'^cm^^ and 10®cm~^, the kinetic temper¬ 
atures between 10 K and 500 K, CO column densities between 
10'®cm“^ to 10^®cm“^ and the beam filling factor rj between 
0.01 and 1. We use the new set of collisional rate coefficients 
calculated by Yang et al. ( 2010) l including energy levels up to J 
= 40 for temperatures ranging from 2 K to 3000 K . A standard 
carbon isotopic ratio '^C/'^C of 70 ( Wilson|[l999 1 is assumed. 
We use the cosmic microwave background radiation at 2.73 K. 

The effect introduced by the default assumption in the cal¬ 
ibration is reflected in a small gap in Eig. between the in¬ 
tensities of the CO lines at the edge of the SSW and SEW re¬ 
ceivers. However this gap is smaller than the 30% error bars as¬ 
sumed for the calculations. Therefore, the inability to use the 
semi-extended calibration tools mentioned in Section lXTI has no 
measurable effect on the temperatures and densities we derive. 


We fit the slope of the cooling curve of the CO lines with 
different combinations of kinetic temperature and gas density, 
which reflects a degeneracy between these two parameters. Eor 
a given combination of kinetic temperature and gas density, we 
obtain the beam-averaged column density by fitting the line ra¬ 


tio *^CO/‘^CO which is sensitive to optical effect^ The optical 
depth at the line center depends on the ratio of the column den¬ 
sity to line width. Eor reasons of simplicity, we take a constant 
width for all the lines measured in the ETS cubes. We assume 
Av = 2 km s“', since for the range and the physical conditions 
considered we find that the RADEX results do not vary signifi¬ 
cantly for Av = 0.5 - 2 km s '. Those values are in agreement 
with HiEi observations towards classical PDRs (e.g., NGC 7023 
Ossenkopf et al.|2013 1 . The beam hlling factor was obtained by 
fitting the mean absolute line intensities. 

The results for three density values, Uh^ = 10^, = 10® 

and nHj = 10®cm“®, are presented in Eigs.[^and|^ and are 
summarized in Table|^ Eor each set of parameters, we derive the 
length of the emission layer along the line of sight: I ~ 

We convert the CO into H 2 column densities taking a standard 


http ://miw. sron.rug.nl/"vdtak/radex/index.shtml 


^ In the cases where '^CO is not detected, we are only able to give a 
lower limit of the column density, assuming a filling factor of one. 
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Fig. 9. Optical depth t determined with RADEX for '^CO (top) and '^CO (bottom) for Hhj = 10"^ cm ^ (left), Hhj = 10^ cm ^ 
(middle) and = lO^cm^^ (right). 


Table 5. The results with RADEX for the different positions. 




nuj 

= 10®cm-3 



HHj = 

10®cm 3 



nH2 = 

10‘'cm 3 



T 

Nco 

1 

7 

T 

Nco 

1 

7 

T 

Nco 

1 

7 


[K] 

(10'® cm- 

-2) (10-3 pc) 


(K) 

(10'® cm- 

-3) (10-3 pc) 


(K) 

(10'® cm- 

-3) (10-3 pc) 


Cl 

70 

100 

4 

0.5 

100 

100 

40 

0.5 

250 

300 

1200 

0.3 

C2 

45 

100 

4 

0.5 

60 

100 

40 

0.5 

100 

200 

800 

0.6 

C6 

40 

> 1 

> 0.05 


70 

> 0.5 

> 0.2 


200 

> 1 

>4 


C8 

33 

> 3 

> 0.1 


42 

> 3 

> 1 


80 

> 5 

> 20 


IRSl 

45 

70 

3 

0.5 

60 

70 

30 

0.5 

90 

200 

800 

0.5 

IRS2 

20 

> 2 

> 0.08 


25 

> 2 

>0.8 


35 

> 5 

> 20 


OFF 

17 

> 3 

> 0.1 


20 

> 3 

> 1 


30 

> 5 

> 20 



relative '^CO abundance to H 2 of 8 x 10 ^ 
let al.|20031 l. 


in PDRs (Johnstone 


Positions Cl, C2, and IRSl are rather similar concerning the 
gas line analysis. We obtain relatively high column densities, 
Aco. of around 1 x 10'* cm“^ of warm and dense gas for Cl, 
C2 and IRSl. The CO temperatures are rather similar in these 
three position while pointing C1 might have slightly higher tem¬ 
peratures. It is rather difficult to fix the density, since with all 
three assumed density values of nn^ = lO"*, Ohj = 10^ and nn^ = 
10® cm~* we are able to fit the observations. The differences lie 
in the obtained length along the line of sight for these three as¬ 
sumptions. Where for nnj = 10"' cm * the length is rather large, 
of around Ipc, for nHj = 10®cm * the length is rather small, 
of around 0.003 to 0.004 pc. The observed projected width and 
extension of the CO emission lines (7 > 9) deduced from the 
detection in the detectors are < 0.1 pc and ~ 0.5 pc, respectively. 
This roughly follows the emission of the dust as seen in Pig. [T] 
The calculated length along the line of sight from the RADEX fit 
is about two times larger than the extension for nn, = 10"' cm“*, 
3 times smaller to the width for nHj = 10® cm“* and 25-30 times 


smaller than the upper width limit for Uhj = 10® cm *. The large 
length along the line of sight for nn, = 10"' cm * exclude lower 
gas densities. The small length along the line of sight derived 
for Ohj = 10® cm~* could be the result of clumps. Regarding the 
comparison of projected extension, 0.1 - 0.5 pc, and obtained 
length from the model calculations, we would therefore tend to 
a density between nn, = 10"' cm“* and nn^ = 10® cm~®. 

Regarding the optical depth, we see from Pig.j^that in point¬ 
ing Cl all the detected '^CO lines are optically thick (E„p ~ 
55 - 500K), while for pointings C2 and IRSl the uppermost 
transitions detected have optical depths t < 1. The '®CO lines 
are mostly optically thin for all positions. 

Por the positions C6, C8, IRS2 and OPE, no '®CO and no 
high-excited '^CO lines are detected, which suggest that the col¬ 
umn density of warm gas is rather small. The obtained lower 
limit on the CO column densities in these four regions are simi¬ 
lar around 1 - 5 x 10'® cm“^. The density and temperature cannot 
be unambiguously determined. 
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Table 6. Line intensity ratios for the central pixel (C3) of each Table 7. Velocity integrated [Cl] (2 - 1)/[CI](1 - 0) ratios and 
pointing. temperatures associated with it for pixel C3 of each pointing. 


Pointing 

CIr 

COr 

CICOtr 

Pointing 

Ratio 

T.y (K) 

Cl 

3.7 ± 1.2 

7.7 ± 0.1 

0.030 ± 0.004 

Cl 

1.6 ±0.5 

138 ± 12 

C2 

3.3 ± 0.9 

4.05 ± 0.08 

0.08 ± 0.01 

C2 

1.6 ±0.4 

131 ± 11 

C6 

1.9 ±0.2 

0.9 ±0.1 

0.41 ± 0.03 

C6 

1.5 ±0.2 

107 ±4 

C8 

1.5 ±0.4 

0.8 ± 0.4 

0.29 ± 0.06 

C8 

1.4 ±0.4 

98 ± 11 

IRSl 

2.4 ± 0.4 

4.64 ± 0.04 

0.08 ± 0.01 

IRSl 

1.5 ±0.3 

114 ±6 

IRS2 

1.3 ±0.2 


1.1 ±0.2 

IRS2 

1.4 ±0.2 

96 ±4 

HII 

1.1 ± 0.2 


2.6 ± 0.5 

HII 

1.4 ±0.2 

91 ±6 

CLay 

1.9 ±0.3 


3.6 ± 0.6 

CLay 

1.5 ±0.2 

109 ±6 

OFF 

1.1 ± 0.2 


0.9 ± 0.2 

OFF 

1.4 ±0.2 

91 ±6 


In position C6 a low density of nnj = 10"*cm“^ results in a 
very high gas temperature. It might be more likely to find lower 
gas temperatures which would lead to slightly larger densities. 

For the positions C8, IRS2 and OFF, it is not possible to 
make constraint on the density. For all assumed densities the gas 
temperatures are smaller than in the other positions. 

For the position OFF we hnd the smallest temperatures, 
which is not that surprising since the position is at a far distance 
from the star and no YSO is assumed in its vicinity. Whether 
this position is located in a dense molecular cloud or in a diffuse 
environment cannot unambiguously be determined. 

For position CLay and HII only the '^CO J=4-3 line is de¬ 
tected, therefore we are not able to carry out a RADEX ht. The 
lack of observed higher-excited '^CO indicates a low tempera¬ 
ture and/or a low density. It can further be assumed that the small 
integrated line intensity of this one line is the result of a rather 
small column density. 

In summary, we get a relatively good assumption of the phys¬ 
ical properties for positions C1, C2 and IRS 1, while for the other 
positions we can only narrow down the physical properties. 


tween regions of RCW 120; group GrA are the YSOs located in 
the PDR, GrB are the dust condensations situated not in the PDR 
but next to it, and in GrC are the pointings targeting regions of 
RCW 120 with diffuse gas. As seen before in section [33] point¬ 
ing IRS2 can be associated with the diffuse gas regions (GrC), 
despite being targeting a known protostellar source. 

In the optically thin limit and assuming LTE, the ratio of the 
velocity integrated emission of [Cl] (2 - 1) and [CI](1 - 0) is 
a sensitive function of the excitation temperature, on the form 
= 38.3K/ln[2.11/C/r] (e.g., [Kramer et al.||2004| l. Using 
this equation, we calculated the temperatures for the central pix¬ 
els of each pointing, shown in column 3 of Table]^ These values 
are different as the ones derived with RADEX. The larger differ¬ 
ences are for the pointings targeting regions off the PDR (IRS2 
and OPE), which lead us to suggest that a large fraction of the 
[Cl] emission is likely originating from the warmer surface lay¬ 
ers. This is also hinted by the [Cl] temperatures, since all of them 
are similar independent of the position in the cloud. 


3.5. PDR model 


3.4. Line ratios 

In Table [^ we give the line ratios CICOtr = 2[CI]/2CO, 
CIr= [CI](2 - 1)/[CI](1 - 0), and COr = CO(8 - 7)/CO(4 - 3) 
for the central pixel of each pointing. In the right panel of Pig.[T0| 
we see how the pointings can be separated into two groups ac¬ 
cording to the COr and CICOtr ratios, labeled GrA (Cl, C2, 
IRSl) and GrB (C6, C8) in the figure. A third group GrC, not 
shown in the figure, can be dehned containing those pointings 
that do not have CO(8 - 7) emission (IRS2, CLay, HII, OPE) 

Por GrA the COr ratio is highest, signaling a high density 
and/or temperature, while the CICOtr ratio is the lowest, sug¬ 
gesting that the cooling by CO lines is more important than by 
[Cl] emission and therefore that the density is high. 

Por GrB, the COr ratio is low and the CICOtr ratio is high. 
The high CICOtr ratio indicates a low density and the low COr 
ratio indicates a low density and/or temperature. 

The pointings in GrC have the lowest COr ratio, since the 
CO(8 -7) line is undetectable, and the highest CICOtr ratio. This 
shows that cooling by [Cl] emission is larger than by CO, and 
much more important than in the other two groups. This medium 
should be therefore less dense compared to GrB. This behavior is 
largely consistent with the temperatures obtained with RADEX 
in Sec. 13.31 

Taking into account the position of each pointing (see Pig. TJ, 
we see that the groups can also be related to a distinction be¬ 


In order to obtain a first approach to the gas density and the 
UV radiation held we applied the Meudon PDR model, devel¬ 


oped by Le Petit et al. (2006 1 . This is a 1-D radiative transfer 
model, which consists of a plane-parallel gas and dust slab of 
a given depth, illuminated on one or both sides by an ultravi¬ 
olet (UV) radiation held and observed in the face-on direction. 
The slab depth is measured by its visual extinction Ay, in mag¬ 
nitudes. The UV held is in units of the local interstellar value of 
5.6 X 10“'“^ ergs cm“^ (Habing held), scaled by a factorPor a 
detailed description of the treatment of the UV radiation held, 
see Appendix C of Le Petit et al. (20061. 

We apply the model to pointing CLay since it is the one that 
better approximates the model assumptions regarding geometry 
and stellar content, and it was set to assume a constant density in 
the slab. We produced models varying the parameters represent¬ 
ing the slab depth (Ay), the atomic Hydrogen initial density (mh), 
and the UV held radiation strength (x), while leaving the remain¬ 


ing parameters in their default settings, as detailed in Le Petit 
[etandSOOhj l. 


We run a grid of models varying log;^' in the range [1,6] 
in steps of 1, lognn from 3 to 6 in steps of 1, and Ay from 
~ 0.5 mag to ~ 500 mag, in tenfold increases. We used the line 
ratios CICOtr = i:[CI]/2COandCIr=[CI](2- iy[CI]a -0)to 
determine the best solution for and;^f. CICOtr is indicative of 
the relative importance of the two main species ([Cl] and CO) 
with regard to the cooling mechanisms. CIr is indicative of the 
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Fig. 10. Line intensity ratios for the central pixel (C3) or each pointing. 
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temperature at the depth in the slab where neutral carbon is lo¬ 
cated and is emitting. 

First, we compare the ratio against the corresponding line 
ratio estimated by the model as a function of hh and;^^, obtaining 
the most likely solution hh ~ 10^; x ~ 10"^- We then test these 
results, running a grid of models this time with x - 10^ fixed, 
and varying tin and Ay. The range of densities is the same as 
before, while Ay ranges between 1 and 10 in unity steps, 
into account the relationship between Nn and Ay (e.g., 
let al.|1978[|Rachford et al.l2002|[EeTetit et al.|2006| ). 


Taking 

Bohlin 


5.34 X 10 


-22 


Nn 


-,-2 


( 1 ) 

( 2 ) 


this means that, by Eq. the width of the slab varies between 
Apdr ~ 6 X 10“^ pc (for log^H = 6; Ay = 1) and Apdr ~ 6pc 
(for log^H = 3; Ay = 10). This is illustrated in Fig. 11 Each 


1.65 X 10^® 

wh 

Apdr 


[cm 

pc 


of the lines represent the locus of points for a given density and 
varying slab width in the [CIr;CICOtr] space. The square is the 
ratio obtained from the data of the central pixel of pointing CLay, 
plotted with its respective error bars. We see that the simulated 
values closest to the observed ratio are the ones for models with 
tin — 10^“^^ and Ay = 4 - 5. 

From the model then, we obtain the values~ 10"^, hh ~ 10^ 
and Apdr between 0.25 and 3.0 pc (with Eq.[^. 


Estimation of ^ at the PDR surface 



log E[CI]/EC0 


Fig. 11. Behavior of tin for different values of Ay, from 1 to 10 
(hlled big square) by steps of 1, for models with;^^ =10"^ plotted 
against the CIr and CICOtr ratios. The line ratios observed at 
each SEW pixel of pointing CLay are plotted with their etTor 
bars. We see that the simulated values closest to the observed 
ratio are the ones for models with tin - 10^ and Ay >4-5. 


from our measurements using the PDR code, especially if we 
consider that dust present in the ionized region, as shown by the 
24 yum emission ( Martins et al.|2010) l, can absorb part of the ra¬ 
diation and diminishes that reaching the PDR surface. 


In order to estimate the UV radiation field impacting on the PDR 
and to compare it with the value derived from the model we 


4. Discussion 


integrated the spectrum of the 08V ionizing star of RCW 120 

(Martins et al. 

2010 

1 between 912 A and 2400 A, following 

Appendix C of 

Le Petit et al. (20061. The star spectrum and 


radius (R* ~ 8.I66R0) were computed and facilitated by F. 
Martins (priv. comm.). The distance of the PDR surface from 
the ionizing star was measured on the SPIRE 350 pn\ image and 
was estimated at ~ 2.15pc. A simple dilution factor (Rtld)^ 
was applied and the backscattering of the radiation by dust at 
the PDR surface was taken into account, as recommended by 
Le Petit et al.| ( 2006) l. We obtain a value of ~ 1925 in Habing 
units. This value is in good agreement with the ~ 1000 derived 


Combining the results from sections 3.3 and 3.5 we see that 
the physical parameters corresponding to RCW 120 are most 
likely between those derived with RADEX with tin - 10"^ - 
10® cm^^, while for the diffuse regions CLay and HII the den¬ 
sity might be lower. Assuming a standard *^CO abundance 
'^CO/Hz ~ 8 X 10 ® iFrerking et Z]|l982| |R611ig et al.||2011 i 
and considering that Ah = A(H) + 2A(H2), we obtain the total 
hydrogen column density values show in columns 3 and 4 of 
Table!] 

On the other hand, using the column density map obtained by 
I Anderson et al.|p012|l for RCW 120, we averaged their column 
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Table 8. A^h obtained for the different pointings. 


Pointing 

Nh 

(xl0^° cm 


Anderson+;2012j 

RADEX (hhj in cm 

10® 10® 10^ 

Nh map 

Cl 

440 

250 

250 

750 

C2 

110 

250 

250 

500 

C6 

100 

> 2.5 

> 1.25 

> 2.5 

C8 

78 

> 7.5 

> 7.5 

> 12.5 

CLay 

12 




Hll 

13 




IRSl 

51 

175 

175 

500 

1RS2 

34 

> 5 

> 5 

> 12.5 

OFF 

28 

>7.5 

>7.5 

> 12.5 


density values obtained at the position of the different pointings 
within the beam of the central pixel of the SLW array (~ 36" 
at 350/im, Makiwa et al. 2013| l, obtaining the total hydrogen 
column densities shown in column 2 of Tabled 

For pointings Cl, C2, and IRSl it was possible to obtain a 
more accurate value for Nh with RADEX and not just a lower 
limit. Among these three positions, only for Cl is the density 
value obtained with Anderson et al. (20121 map between the 
RADEX values. Nevertheless, the column densities derived by 
both methods are, in a first order, in agreement. The assump¬ 
tions on the dust properties play a crucial role when deriving 
the column density from dust observations, therefore we cannot 
make a detailed comparison. 

Eor the other pointings, with RADEX we only obtain lower 
limits to the column density, and correspondingly, the average 
values from Anderson et al. ( |2012| l map are higher than the 
RADEX values. In particular, pointing C6 shows a large differ¬ 
ence between the two estimations, in line with the suggestion 
made in section 3.3 that its temperature is lower than that ob¬ 
tained with RADEX. 

Any of the two methods show the differences in gas proper¬ 
ties throughout RCW 120. The similar values for average col¬ 
umn density of pointings CLay and HII suggest that the PDR of 
RCW 120 also extends on the plane-of-the-sky direction, sup¬ 
porting the bubble-shaped morphology proposed for it. 

Pointings C6 and C8 show temperatures and densities inter¬ 
mediate between the warm and dense regions on the PDR, and 
its diffuse and cooler parts. As mentioned in section [33] the rel¬ 
atively high temperature given by the code for pointing C6 is 
likely overestimated, leading to a slightly larger densities. 


5. Summary 

We have obtained Herschel SPIRE-ETS spectra towards 9 po¬ 
sitions in the RCW 120 Hii region, detecting the [Cl] lines at 
370 and 609 yum, the 205/rm [Nil] transition, the ^^CO ladder 
between the 7 = 4 and 7=13 levels and the '^CO ladder be¬ 
tween the 7 = 5 and 7 = 14 levels. CH^ was detected in ab¬ 
sorption at all positions in the region, however the low spectral 
resolution of the spectra do not allow us to obtain quantifiable in¬ 
formation from that line. Nevertheless, its ubiquitous absorption 
suggest the presence of diffuse gas along the line of sight, while 
RCW 120 may emit in this line with varying absorption depth. 
The [Nil] emission line is strong and is also detected over the en¬ 
tire field. The [Cl] lines are detected in almost all detectors with 
a ratio which shows little variation throughout the region. This 


suggests the presence of low-density PDR over the entire RCW 
120 region. This is further supported by the temperatures ob¬ 
tained with the ratio of the two [Cl] lines detected, which show 
little variation throughout the region. 

The low-excitation '^CO lines are detected everywhere, 
while higher-excited lines are only detected in the condensations 
Cl, C2 and IRSl together with '^CO. We use RADEX to derive 
the physical properties at these positions. The gas temperatures 
are 45 - 250 K for densities of 10"^ -10® cm“^, and a high column 
density that is in agreement with dust analysis. The excited CO 
could arise either from the edge of the dense irradiated structure 
or small dense clumps containing young stellar objects. We see 
the excited CO emission in several detectors, partly in an elon¬ 
gated region, coming from the PDR and/or several young stellar 
objects. The analysis of the other condensations C6 and C8 re¬ 
veal a less dense medium with still high gas temperatures. Eor 
the positions HII, CLay and OEE, where no condensations are 
observed, reveal the lowest densities with a highest CICOtr ra¬ 
tio. 

We model the PDR of RCW 120 (poiting CLay) with the 
Meudon PDR code. We obtain a hydrogen density of hh ~ 
10^'^cm^^ and an ionizing radiation fields ~ 10^ in Habing 
units. The value for x agrees with what is expected from the 
emission of an 08V star, as is the ionizing star of RCW 120. 
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